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Using 58 million J/xp events collected with the Beijing Spectrometer (BESII) at the Beijing 
Electron-Positron Collider, the decays J /tp "/(pp and J/ip yup are searched for, and upper 
limits on their branching fractions are reported at the 90% C. L. No clear structures are observed in 
the 7p, 70, or ptp mass spectra for J/tp — > ytpp nor in the 7p, yui, or pcj mass spectra for J/ip — > fup. 
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PACS numbers: 13.20.Gd, 13.25.Gv, 13.20.-v, 12.38.Qk, 14.40.-n 



I. INTRODUCTION 

QCD predicts a rich spectrum of gg glueballs, ggq hy- 
brids and qqqq four quark states along with the ordinary 
qq mesons in the 1.0 to 2.5 GeV/c^ mass region. Radia- 
tive J/ip decays provide an excellent laboratory to search 
for these states. Until now, no unique experimental sig- 
natures of such states have been found. 

Systems of two vector particles have been intensively 
studied for signatures of gluonic bound states. Pseu- 
doscalar (0~) enhancements in pp, uiui, and (jxj) final states 
have been seen in radiative J /ip decays [1-7] , and a scalar 
(0+) enhancement near ucj) threshold is observed from the 
doubly OZI suppressed decay of J/V' ^uxj) with mass 
M = 18121^^ ± 18 MeV/c^ and width T = 105 ± 20 ± 28 
MeV/c^ [8]. The radiative J /if) decays J/tj) ^p<j) and 
J/tp — > are OZI suppressed processes, and the mea- 
surements of these two decays and the search for possible 
resonant states in their decay products will provide useful 
information on two vector meson systems. 

The double radiative channels J/xj) ^ 'yX, X jV 
{V=p, (f) and w) are studied to probe the quark content 
of the object X. The r/(1440) has been studied by the 
BES Collaboration through the double radiative channels 
J/tp ^ jiip) and J/ip ^ lilfp) [9]. At one time, the 
Ty(1440), after it was observed in J/i/; decay [10], was 
regarded as a glueball candidate. But this viewpoint 
changed when its radiative decay modes were observed. 
The 7?(1440) is seen (at 1424 MeV/c^) by the BES Col- 
laboration to decay strongly into 7p, not ^(f). Prom this 
result, one cannot draw a definite conclusion on whether 
the r;(1440) is either a qq state or a glueball state. There- 
fore, further study is needed to clariiy the situation. The 
process J/i]j VX, X — > -yV also allows us to study X 
using the system. 

In this letter, we report on the measurements oiJ/ip^ 
"fp(f> and J/ip — > "fpuj decays and the search for possi- 
ble structure in the "fV and VV invariant mass spectra, 
using 58M J/ip events collected with the Beijing Spec- 
trometer (BESII) at the Beijing Electron-Position Col- 
lider (BEPC) [11]. 



II. DETECTOR AND DATA ANALYSIS 

BESII is a conventional solcnoidal magnet detector 
that is described in detail in Refs. [11]. A 12-layer vertex 
chamber (VC) surrounding the beam pipe provides coor- 
dinate and trigger information. A forty-layer main drift 
chamber (MDC), located radially outside the VC, pro- 
vides trajectory and energy loss (dE/dx) information for 
tracks over 85% of the total solid angle. The momentum 
resolution is ap/p = 0.017-\/l + (p in GeV/c), and the 
dE/dx resolution for hadron tracks is ^8%. An array of 
48 scintillation counters surrounding the MDC measures 



the time-of-flight (TOF) of tracks with a resolution of 
'~200 ps for hadrons. Radially outside the TOF system 
is a 12 radiation length, lead-gas barrel shower counter 
(BSC) which measures the energies of electrons and pho- 
tons over ~80% of the total solid angle with an energy 
resolution of aE/E =22%//E {E in GeV). Outside of 
the solenoidal coil, which provides a 0.4 tesla magnetic 
field over the tracking volume, is an iron flux return that 
is instrumented with three double layers of counters that 
identify muons of momentum greater than 0.5 GeV/c. 

A. General selection criteria 

All the charged tracks are reconstructed in the MDC, 
and the number of charged tracks is required to be four 
with net charge zero. Each track should (1) have a good 
track fit; (2) have | cos6\ < 0.8, where is the polar an- 
gle of the track measured by the MDC; (3) originate from 

the interaction region, ■^V^ + Vy < 2 cm and \ Vz\ < 20 

cm, where Vx, Vy, and Vz are the x, y, and z coordi- 
nates of the point of closest approach of the track to the 
beam axis; and (4) be identified as either a pion or a 
kaon. The particle identification (PID) is accomplished 
using the TOF information and the dE/dx information 
from the MDC [9]. For instance, a pion should have a 
higher particle PID confidence level than those for other 
hypotheses (kaon, proton). 

A neutral cluster in the BSC is considered to be a 
photon candidate when its energy deposit in the BSC 
is greater than 30 MeV, the angle between the nearest 
charged track and the cluster is greater than 10° (to reject 
photons associated with the charged particles), the first 
hit appears in the first five layers of the BSC (about 
six radiation lengths of material), and the angle between 
the cluster development direction in the BSC and the 
photon emission direction is less than 25°. If the angle 
between two photons is less than 15°, and the invariant 
mass of these two photons is less than 60 MeV/c^, the 
two photons are merged. 

B. Analysis of J/ ip ^ l4>P 

For the channel J/'ip — > {(j) — > K^K~, p — > 
TT+TT"), we require two pions, two kaons, and at least 
one photon with energy greater than 50 MeV. Next, the 
selected charged tracks and photon are fitted kinemat- 
ically using energy and momentum conservation con- 
straints (4C), looping over all photon candidates. The 
combination with the minimum is selected, and the 
photon in this combination is taken as the radiative pho- 
ton. Next, we require the energy of the radiative photon 
to be greater than 0.1 GeV and the fit to be less than 
8.0. The requirement is determined by optimizing the 
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signal (5) to noise (B) ratio {S/\/STB), where S + B\s 
determined from data and S is determined from Monte 
Carlo (MC) simulation after event selection. 

In order to suppress backgrounds from J/ijj two body 
decay channels with a ^, for example J/tp (f)rj and 
J/^/i — > (^77' (958), the momentum is required to satisfy 
P4, < 1.1 GeV/c. To select (j> particles, the K+K' in- 
variant mass is required to satisfy \Mx+k- ^ 0-01 
GeV/c2. 

The main expected background channels can be di- 
vided into the following five groups: (1) J/ij} — > (pr], 
(l)V'{958); (2) J/tP uoKK, <f>KK -> tt+tt-tt"); 
(3) J/^ ^ a.-/o(980) (/o(980) ^ K+R-); {^)_JJ^ ^ 
0./()(98O): and (5) J/i) {-f)K*{892)K*{892), 
K^K~ . Since J/ij} ^ TT^p(j) is forbidden by 
C-parity conservation, it can be neglected in our back- 
ground analysis. The tt+tt" invariant mass distributions 
from all above possible background channels are smooth 
with no p peak according to MC simulations; therefore 
they will not affect the determination of the number of 
signal events. 



(7, 7r'')7r^7r 





2: (a) The yp invariant mass spectrum after (p selection 
M^\ < 0.01 GeV/c^), (b) the j(f) invariant mass 



1.0 1.2 1.4 

M^+,- (GeV/c=) 

FIG. 1: The tt^tt" invariant mass distribution. The square 
points with error bars are data, the shaded histogram is from 
the <j) sideband events, the dotted curve is the signal shape 
from MC simulation, and the blank histogram is the fit. 

The TT+TT^ invariant mass distribution for events that 
survive the selection criteria is shown in Fig. 1, and 
the TT+TT" invariant mass distribution from (p sidebands 
(1.05 < Mk+k- < 1-08 GeV/c2 or 0.985 < Mk+k- < 
0.992 GeV/c^) scaled to the amount of background in the 
signal region is also shown as the shaded histogram. 

By fitting the tt+tt" invariant mass distribution with 
a p signal shape obtained from MC simulation and using 
the histogram from cj) sideband events to describe the 
background shape, 43.2 ± 18.8 signal events are obtained, 
as shown in Fig. 1. The p signal statistical significance 
is estimated by comparing the likelihood values with and 
without the signal in the fit, and it is only about 2a. The 
detection efficiency is (3.33±0.04)% from MC simulation. 

Figures 2 (a) and (b) show the invariant mass distri- 
butions of 7p and ^(p; and Figs. 2 (c) and (d) show the 



FIG 

{\Mk+k- .... 

spectrum after p selection {\M^+^- — Mp\ < 0.15 GeV/c^), 
(c) the pcj) invariant mass spectrum, and (d) the Dalitz plot 
of M^^ versus M^p. 



pcj) mass distribution and Dalitz plot of M^^ 

for J/tp jcpp candidates, where, the tt'^tt' 
mass of the p candidates must satisfy |M^+^- 



versus M^^ 
invariant 
-MJ < 



0.15GeV/(r. No clear structure around 1440 MeV/c^ 
region is observed in the M^p distribution. There may 
be some structures in the and p4> mass spectra, but 
because of the low statistics, it is difficult to determine 
whether they are real resonances or just statistical fluc- 
tuations. 

The systematic errors, are evaluated with selected sam- 
ples that are compared with MC simulations. In this 
analysis, the systematic errors on the branching fraction 
mainly come from the following sources: 

(1) . Particle identification (PID). In Ref. [12], the 
PID efficiencies of pions and kaons are analyzed in detail. 
Here, two charged tracks are required to be identified as 
pions, and the other two are required to be kaons, so the 
systematic error from PID should be less than 4%. 

(2) . MDC tracking and kinematic fit. In order to study 
the systematic errors from the MDC tracking and kine- 
matic fit, many distributions from data, including the 
wire efficiency and space resolution of hits in the MDC, 
are compared with those from MC simulations, using two 
different treatments of the wire resolution simulation. 
The difference between the two simulations is taken as 
the systematic error for the tracking [12]. In this paper, 
13.5% is conservatively taken to be the systematic error 
from MDC tracking and kinematic fit. 

(3) . Photon detection efficiency. The photon detection 
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efficiency is studied using J/ip p^n^ in Ref. [12]. The 
results indicate that the systematic error is less than 2% 
for each photon. 

(4) . Background uncertainty. In order to determine 
the background uncertainty, several background assump- 
tions were tried, including: (a) a fourth order polynomial 
function, (b) the background shape from MC simulations, 
and (c) the histogram from (j) sideband events. The dif- 
ferences between the different fit results are taken as the 
systematic error, which is about 13.0%. 

(5) . Intermediate decay branching fraction and the un- 
certainty of the number of J /ip events. The <f) decay 
branching fraction (1.2%) from Ref. [13], and the un- 
certainty in the total number of J/V' events (4.72%) are 
also considered as sources of the systematic error. 

Adding all errors in quadrature, the total error is about 
20.0%. 

Finally the branching fraction is: 

Br{Jltl) 7#) = (4.5 ± 2.0 ± 0.9) x 10-^ 

where the first error is statistical and the second is the 
systematic. Since the statistical significance of the p sig- 
nal is only 2a, the upper limit (90% C.L.) is also deter- 
mined by a Bayesian method [13]: 

Br{J/tp ^(pp) < 8.8 X 10"^ 



C. Analysis of J/V" — > 7wp 



7r"'"7r tt^jP 



For the channel J/tl) ^ujp {ui 
7r+7r^,7r" 77), we require four pious and greater 
than two photons, where the energy of the photon candi- 
dates should be greater than 50 MeV. Next, the selected 
charged tracks and three photons are fitted using a 4C 
kinematic fit, looping over all photon candidates, under 
the hypothesis of J/ij} — > 372(7r+7r~). The combination 
with the minimum is selected, and the is required 
to be less than 9.0, which corresponds to the best signal- 
noise ratio. The two photons with invariant mass closest 
to the 7r° mass are regarded as being from 7r° decay, and 
the other is taken as the radiative photon. Finally, a 
5C kinematic fit is made under the J/-^ ^ 72(71 "'"7r~)7r° 
hypothesis with the invariant mass of the 77 pair asso- 
ciated with the 7r° being constrained to m^o. After the 
5C kinematic fit, we require to bo less than 9.0 and 
the energy of the radiative photon to be greater than 0.1 
GeV. 

The 7r"'"7r~ and 7r~*'7r~7r*^ combinations for the p and UJ 
candidates are selected from the minimum value of 



where Cp and arc the widths (the mass resolutions are 
included) of p and ui, respectively, from MC simulations. 



To select the uj candidates, we require \M^+^-^o—M^\ < 
0.023 GeV/c2. 

In order to suppress the backgrounds from J/ip two 
body decay channels with an lu, for example J/iJj ujrj 
and J/if) LOT]' {958), the w momentum must satisfy 
Pa; < 1-2 GeV/c. Similarly we also use Pp < 1.1 GeV/c 
to suppress the backgrounds from J/V' two body decay 
channels with a p resonance. 

The possible background channels of J/ijj "'MP can 
be divided into the following four groups: (1) J/ip — > 
ujr], ojri'{958), wtt+tt^, 6^/0(980), -yojuj, 6i(1235)±7r=F; 

(2) J/ip a2(1320)p, pr], pi, -/pp, 777^ [Vc ^ Pp); 

(3) ,J/ip — > 2(7r+7r^)7r'', 72(7r+7r~), 'yrjTnr, ^r(; and 

(4) J/V' ^ 7X (X -> ?77r7r, K*'W). Since J/V' ^ 77° 

is forbidden by C-parity conservation, it can be neglected 
in our background analysis. The contamination from all 
the above possible backgrounds after event selection is 
found to be very small according to MC simulation. 

Figures 3 and 4 show tt^tt^tt'^ and tt+tt^ invariant 
mass distributions and the tt+vt^ mass distribution of lo 
sidebands (0.6 < M^+^-^o < O.STGeV/c^ or 0.915 < 
Mjr+^-^o < 0.95GeV/c^), scaled to the amount of back- 
ground in the signal region and shown as the shaded 
histogram in Fig. 4. After subtracting the w sideband 
events, we fit the tt+tt" invariant mass distribution to 
obtain the number of J/iJj ^ojp events. 

Although the branching fraction of p ^ tt+tt" 
(-^100%) is about two orders of magnitude larger than 
that of a; — > tt+tt" (1.7%), the interference between p 
and ui must be considered, and therefore, the fit function 
is expressed as [14]: 



N{M„+^-) = L{M^+^-)+\Ap{M^+^-)+A^{M^+^-)e 



where i is a polynomial background term, </? is the rel- 
ative phase angle between the two amplitudes, Ap and 
A^, which are represented by Breit-Wigner functions up 
to a numerical factor: 

Av{M^+^~) = y^Fewvi^^+^-W = P^)^ 

and the Breit-Wigner functions have an s-independent 
width: 



Fbw 



TM 



s-M'^+iMT' 



Finally, the number of p events (181.3±72.7), the num- 
ber of w events (76.6 ±61.2), and the relative phase angle 
(39.4±36.4)° are obtained; the fit is shown in Fig. 5. The 
statistical significance of the p signal is about 3.1cr, and 
the detection efficiency is about (1.36 ±0.01)% from MC 
simulation. 

Figures 6 (a) and (b) show the invariant mass distribu- 
tions of ju) and jp, and Figs. 6 (c) and (d) show the pui 
mass distribution and Dalitz plot of versus M^^ for 
the J/tp ^iop candidates, where the p signal is selected 
using |M^+^- - Mpl < Q.lbGeV/c?. No clear structure 
around 1440 MeV/c^ region is observed in the M^p dis- 
tribution. Although there is a hint of a possible structure 
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FIG. 3: ThoTT 7r~7r invariant mass spectrum. 
The uj signal and u) side band regions are in- 
dicated by arrows. 




0,4 0,6 0,8 1 1,2 1,4 1.6 
M.+.- (GeV/c^) 

FIG. 4: The ■k'^ti~ invariant mass distribu- 
tion. The shaded histogram is the distribution 
from u) sideband events which has been scaled 
to the amount of background in the signal re- 
gion. 



around 1700 MeV/c^ in the pu mass spectrum compared 
with the scaled puj mass distribution from p and u side- 
band events, it is difficult, because of the low statistics, 
to determine whether it is a real resonance or just due 
to a statistical fluctuation. Finally, only the branching 
fraction of J/ip jivp is given. 

Systematic errors in the J/tp ^ jLOp branching frac- 
tion measurement arc analyzed similarly as in the J/ip 
^(pp channel, which mainly come from particle identifi- 
cation (4%), the MDC tracking and the kinematic fit 
(11.5%), the photon detection efficiency (6%), the ffiting 
procedures and different treatments of the backgrounds 




M^+^- (GeV/c" 



FIG. 5: The fit, described in the text, to the Tv'^n^ invariant 
mass distribution. Here the squares with error bars arc data, 
the dotted curve is the total signal from p and uj (the inter- 
ference between them is also included), the dashed curve is 
the polynomial background, and the solid curve is the fit. 
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FIG. 6: (a) The 70; invariant mass spectrum after p selec- 
tion i\M^+„- - Mp\ < 0.15 GeV/c^), (b) the 7p invariant 
mass spectrum after u) selection (|M^+^-^o — M^^] < 0.023 
GeV/c^), (c) the puj invariant mass spectrum (the shaded his- 
togram is the scaled pu; mass spectrum of p and ui sidebands 



events), and (d) the Dalitz plot of M^p 



versus M^^. 



(21.6%), the total number of J/^ events (4.72%), and the 
UJ decay branching fraction (0.8%), taken from Ref. [13]. 
The total systematic error is 26.0%. 

Finally, the branching fraction can be obtained: 

Br{ J/iP jLup) = (2.6 ± 1.1 ± 0.7) x 10"^, 

where the ffist error is statistical and the second is the 
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systematic. Since the statistical significance of the p sig- 
nal is only 3.1cr, the upper hmit (90% C.L.) is also esti- 
mated by a Bayesian method [13]: 

Br{J/ij jLup) < 5.4 X 10"^. 

III. CONCLUSION 

Table I summarizes our results for the J/tp jVV 
branching fractions; it also lists our J/tp ^ Jujw branch- 
ing fraction with one of lu decays to tt~^tt~ , which has a 
very small branching fraction (15.9%). Taking into con- 
sideration its large statistical error, this branching frac- 
tion is consistent with the value from Ref. [13]. 

As a check, the J/ip jujuj branching fraction was 
fixed at the value from Ref. [13] and the previous result 
from BES Collaboration [1], and the fit was redone to ob- 
tain the J/tp ^ jcop branching fraction. This branching 
fraction is not sensitive to that of J/iJj ^lulo. 

From our analysis of these two channels, we did not 



observe any clear structures in either the VV or mass 
spectra. 
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TABLE I: The branching fractions of J/i) ^VV 



Doca. 


V Alodo 






BESII' 




— > 'yulp 


(2.6 ± 1.1 ± 0.7) X 10"" 
(< 5.4 X lO-"") (90% C.L.) 






J/i> 


i4>p 


(4.5 ±2.0 ±0.9) X 10"" 
(< 8.8 X 10-5) (90% C.L.) 






J/V 
one uj 


— > 7r"'"7r~ 


(6.0 ± 4.8 ± 1.8) X 10"^ 
(< 1.7 X 10-2) (90% C.L.) 


(1.59 ± 0.33) X 10-^ 


(2.29 ± 0.08) X 10-^ 



t Result from Jjii ->■ 777(1770) ->■ jujuj [1] 



